Charge oscillation-induced light transmission through subwavelength slits and holes 
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We present a concrete picture of spoof surface plasmons (SSPs) combined with cavity resonance to clarify the 
basic mechanism underlying extraordinary light transmission through metal films with subwavelength slits or 
holes. This picture may indicate a general mechanism of metallic nanostructure optics: When light is incident 
on a non-planar conducting surface, the free electrons cannot move homogeneously in response to the incident 
electric field, i.e., their movement can be impeded at the rough parts, forming inhomogeneous charge distribu- 
tions. The oscillating charges/dipoles then emit photons (similar to Thomson scattering of x rays by oscillating 
electrons), and the interference between the photons may give rise to anomalous transmission, reflection or 
scattering. 

PACS numbers: 42.25.Bs, 42.79.Dj, 78.67.-n, 84.40.Az 
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Since the discovery of extraordinary light transmission 
through metal films perforated by subwavelength hole arrays 
Ql], tremendous theoretical and experimental work has been 
carried out to understand the underlying physics. Several 
mechanisms, particularly surface plasmons (SPs), have been 
proposed as the possible origins S 0, HI S 13, H] • How- 
ever, no universal understanding has been reached to date. 
Here based on first-principle calculations, we present a charge 
oscillation-induced light emission mechanism, which gives 
the origin of enhanced transmission in subwavelength systems 
and may also shed light on the fundamental of interactions be- 
tween light and metallic nanostructures. 

In non-magnetic media, the electric and magnetic fields are 
coupled by Maxwell's equations V x E = —iKH and V x 
H = iKeEi (K = 2tt/X, A the wavelength in vacuum, and 
e the effective permittivity). For varying e(r) (i.e., Ve ^ 0), 
the divergence of the second equation generally gives V ■ E = 
— [(Ve) • E]/e = 4irp ^ 0, where p is the charge density. 
Consider in Fig. 1 the free-standing one-dimensional (ID) 
gold grating illuminated by a plane wave, where Maxwell's 
equations can be solved by the rigorous coupled-wave analy- 
sis (RCWA), a first-principle method |9[]. For simplicity, we 
only discuss normal incidence in this Letter. Figure 1 shows 
the zero-order transmittance (Tq) spectra of both TM (H || y) 
H and TE (E || y) waves fl(J]- Here RCWA correctly re- 
veals the cutoff wavelength A c ~ 2W, above which transmis- 
sion of TE waves is forbidden. The reason is that TE waves 
in the slit approximately take the waveguide modes E y oc 
sm{miTx/W) exp[± mz{A/\ 2 — m 2 /W 2 ) 1 ^ 2 ] (m/0 being 
integers) flTll . where for A > 2W, all the modes are evanes- 
cent. The drastic difference between the TE and TM spectra 
stems from the fact that TE waves satisfy V • E = while TM 
waves may induce electric charge oscillation p(r)e lut (u> the 
frequency), where p(r) = V • E(r)/47r. 

For TM polarization, the vertical component E z (x, z) of 
the electric field (invariant with y) has abrupt discontinuity 
across the surfaces at z s = and h, from which one obtains 
the surface charge density p s (x,z s ) = SE z (x, z s )/47r. In 



Fig. 2(a), the calculated p s (x, 0) profile correctly shows that 
charges only exist on the metal surface [p s (x, 0) = for < 
x < W], and the charges tend to accumulate at the metal cor- 
ners. At resonant wavelengths, p s (x, 0) ~ (— l) N p s (x, h), 
i.e., the charge patterns on the two surfaces are nearly the 
same, but they have opposite signs for odd resonant num- 
bers N (defined below). The bulk charge density is given by 
p v = V ■ E/4-7T. In Fig. 2(b), the p v (x,z = const) curve 
calculated from the internal E eigenmodes reveals that there 
is a strong peak exactly centered at each slit wall, x w = 
or W (plus any multiple of d). When sufficient diffrac- 
tion orders (1601 orders in Fig. 2) are retained in RCWA, 
p v (x, z = const) approaches a delta function across x w [see 
Inset of Fig. 2(b)], i.e. charges inside the grating also appear 
as surface charges on the walls. 

The profile p v (0,z) [= —p v (W,z)] plotted in Fig. 2(c) 
shows that in the central range the charge density on the wall 
is nearly a standing wave (with wavevector k z ~ 2w/X). 
This indicates that the electric field in the slit is also a 
standing wave lfl2l [l3ll consisting of a forward wave E a ~ 
E a exp(—ik z z)x and a backward wave E& ~ Eb exp(i/c z z)x 
in Fig. 3. Note that \p v (x w , z) \ increases sharply when z — > z s 
in Figs. 2(c). Mathematically, the charge density at each cor- 
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FIG. 1: (color online) Transmission spectra of a gold grating with 
period d — 3.5, slit width W = 0.5, and thickness h = 4 pm. x, 
y, z are unit vectors. The TM' curve was calculated with Re[£ c (u>)] 
(< 0) for gold replaced by — Re[e c (a>)], 
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FIG. 2: (Color online) Charge distributions in the gold grating of 
Fig. 1 under resonance. TM polarization. A = 4.845 pm (for peak 
N — 2 in Fig. 1). (a) Surface charge densities, (b) Sectional charge 
densities, (c) Charge densities on the wall. 



ner consists of both p s (x w ,z s ) and p v (x w ,z — ► z s ). Since 
they are always in phase, these two contributions are superim- 
posed constructively. This makes the total charge densities at 
the corners much higher than in other regions, which is a typi- 
cal edge effect. Consequently, there exist two large oscillating 
dipoles at the two ends of the slit. 

Now we may obtain a clear picture about the light scatter- 
ing process. As shown in Fig. 3, the incident electric field Ej„ 
drives (mainly) free electrons on the metal surface to move, 
but the movement is blocked at one corner of the metal, re- 
sulting in accumulation of electrons there. Meanwhile, ex- 
tra positive charges appear at the other corner since some of 
the electrons have moved away. Thus, two dipoles P a and 
P r are formed. Since they oscillate with the incident wave, 
these dipoles act as light sources emitting new wavelets (pho- 
tons), which form scattered waves (Thomson scattering). If 
we consider each period of the array as an overall scatter- 
ing unit, the wavelets emitted from two adjacent units have 
a path difference d sin 9 along an arbitrary direction 9. Note 
that for a subwavelength slit array with d < A, we have 
dsm9 < d < A, so the phase difference can never reach 
2ir. This means that the oblique wavelets always tend to can- 
cel each other out in the far fields (destructive interference, 
similar to the absence of x-ray diffraction at non-Bragg an- 
gles). Thus, they form evanescent waves near the surface. 
Mathematically, these evanescent waves can be expressed as 
E m exp[-(G^ - K 2 y/ 2 \z\ - iG m x], where G m = 2rmr/d 
(m 7^ being integers) and K = 2tt/\ (note that \G m \ > K 
for d < A) in. 

One may prove that for any wavelength A > d, the charge 
patterns on the upper surface are always the same as that in 
Fig. 2(a) except that the peak heights vary with A. There- 
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FIG. 3: Light transmission process in the ID grating. The phases 
(directions) of the waves and dipoles are based on Fig. 2 (and they 
all have a common oscillating factor e ). 



fore, the period of the charge density wave always equals 
the lattice constant d and generally does not satisfy the dis- 
persion relation k sp = K[e c /(l + Ec)] 1 ^ 2 of a classical SP 
(CSP) [ 14] at all (e c the conductor's permittivity and k sp the 
wavevector of the CSP). However, the evanescent waves are 
indeed surface-bound modes with their strengths decaying ex- 
ponentially along —z. So they are spoof SPs (SSPs) Q1, 
but their formation is the result of charge oscillation-induced 
light emission and destructive interference. The SSP model 
has been conceptually proposed by Pendry et cd. Q3[] (also see 
fi^l ). Here we provide a concrete picture illustrating its ori- 
gin. 

The wavelets propagating along the backward direction 
(9 = 0) are always in phase (the path difference is zero). 
Therefore, the back reflected wave is not evanescent but a 
propagating mode. However, since E # and E r always have 
opposite directions, they tend to offset each other, thus reduc- 
ing the overall back reflection. 

In Fig. 3, dipole P a also emits a wavelet E a inside the slit. 
Similarly, electrons on the slit walls also move in response 
to the electric field in the slit (represented by the current J). 
Then the charge movement is disrupted again at the exit edges, 
giving rise to another large dipole P^. Through this tunneling 
process, the charge patterns on the upper surface are dupli- 
cated on the lower surface. Oblique wavelets emitted from 
these duplicated light sources then form SSPs again below 
the grating. The sub-wavelets Et and E t along the forward 
direction form the zero-order transmitted wave. Overall, the 
grating only emits two propagating modes, the reflected and 
transmitted waves, that share the incident energy, while all the 
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other modes are SSPs. 

The oscillating dipole can give a strong feedback to the 
upper surface by emitting a wavelet propagating upward. 
If Eb is not in phase with E a (and Ej„) at z = 0, it sup- 
presses the strengths of P a and Er. Then E r , which includes 
specular reflection from the metal surface, becomes domi- 
nant, leading to strong backward reflection. The weakened 
E meanwhile reduces P;,. However, if E& is in phase with 
E a at z = 0, it enhances P a . The enhanced P a subsequently 
strengthens E a , P&, E&, and so on. Then a Fabry-Perot-like 
resonant state is formed with the strengths of all the dipoles 
and wavelets maximized. Under this condition, E r is largely 
offset by E^ in the far-field regions, leading to minimized 
backward reflection. Below the grating, wavelet E t also partly 
offsets Et, but |Ex| can be much larger than |E t | (unlike 
E r , E t does not include any specular reflection). Therefore, 
when Et is maximized, it maximizes To. The resonant wave- 
length is always slightly longer than 2h/N, the ideal Fabry- 
Perot wavelength, where N is the resonant number (number 
of the standing wave nodes). Although the Fabry-Perot-like 
resonance has been recog nized before fl & H M d H, 
here we explictily illustrated its origin. Particularly, Figure 3 
shows that it is dipoles P Q and P;, that act as the two "reflect- 
ing planes" required to form a vertical resonator 1 1311 . 

Note that when A < d, some of the oblique waves satisfying 
\Gm\ < K become propagating diffracted waves (i.e., they 
are no longer SSP modes). The sharing of the incident energy 
by these diffracted waves significantly reduces To (and back- 
ward reflection). This explains the much lower transmittance 
in the A < d range in Fig. 1. 

Light scattering from a 2D hole array has a similar pic- 
ture. The incident electric field drives electrons on the upper 
surface to oscillate. The charge movement is blocked at the 
hole edges, giving rise to oscillating dipoles at the entrance 
openings. Thus, the holes act as a 2D array of light sources. 
For a subwavelength lattice with max(<i, c^) < A (Fig. 4), 
the phase difference between wavelets emitting from adja- 
cent holes again is less than 2tt along any oblique directions. 
So these wavelets form SSPs above the film except that the 
wavelets along — z constitute a non-evanescent reflected wave 
(which offsets specular backward reflection). Meanwhile, the 
charge patterns are tunneled onto the lower surface due to the 
charge movement on the hole walls. Consequently, a similar 
set of SSPs and a forward transmitted wave are formed below 
the film. 

However, 2D holes have a different tunneling mecha- 
nism. Consider the rectangular hole (a unit cell of an ar- 
ray) in Fig. 4. Except for the wave distortions near the 
ends, the electric fields in the hole roughly take the waveg- 
uide modes. For subwavelength holes, the basic 2D TE^o 
modes dominate JH [01]. When A > 2L, the modes in the 
hole are evanescent and approximately take the forms E a oc 
sm(iry/L) exp(— /3z)x and Ef, oc suijiry/L) exp|— (3(h — 
z)]x, where (3 = n(l/L 2 - 4/A 2 ) 1/2 |Qj]. E drives charges 
on the walls to move/oscillate, but the charge density de- 
cays with decaying E a toward +z due to the waveguide re- 
striction. Nevertheless, the sharp discontinuity of the charge 
movements at the exit end can still cause significant accumu- 




FIG. 4: Light transmission through 2D hole arrays, (a) The unit 
cell of 2D lattice viewed along z. (b) Side view. Here the dipoles 
similar to those denoted by P r in Fig. 3 are omitted, (c) Charge 
densities \p(x, y, 0) | of a silver hole array: h = 100, W = L = 150, 
d — di = 600 nm. A = 673 nm, corresponding to peak P in the 
zero-order transmittance (Too) spectrum in Inset I. Inset II: Profiles 
p(0, y, 0) and p(0, L/2, z) with constant phases at A = 673 nm. 



lation of charges there, giving rise to a relatively large dipole 
Pb, provided that h is adequately small. Since E Q and E^ 
no longer have position-dependent phase factors, they are al- 
ways in phase with each other and the incident wave, and 
they always resonate as the feedback E^ is always construc- 
tive. But the strengths of the dipoles and wavelets decrease 
with increasing h due to the decaying effect. For fixed h 
and in the absence of diffraction effects, the transmittance 
would increase monotonically with A decreasing toward 2L 
(f3 decreasing). However, when A < maxfrf,^), diffracted 
waves emerge, which reduces the zero-order transmittance in 
the short- wavelength range. Therefore, strongest transmission 
must occur for A > max(c?, c^), the non-diffraction range, but 
meanwhile A should still be close to max(d, c^) so that the 
damping coefficient (3 is sufficiently small. 

This picture agrees very well with the measured transmis- 
sion spectra in the literature, and our finite-difference time- 
domain (FDTD) computations have unambiguously proved 
it. Figure 4(c) shows the calculated charge distribution 
\p(x,y,0)\ on a free-standing silver film with a hole array, 
where the charges again exist mainly on the two walls x = x w 
(perpendicular to the charge movement direction). The resid- 
ual charges in other regions are surface charges and they 
disappear for < z < h. (Here note that FDTD gives 
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the overall charge densities.) The shape of the p(0, y, z) 
[= —p(W,y,z) ~ sin(7ry/L)] profile is independent of z, 
but the maximum density p(0, L/2,z) changes with z in inset 
II. Here our calculations indeed reveal the two large dipoles at 
the two openings of the hole (also see IU9I1 for the single-hole 
case). 

The above illustrations may indicate a general Thomson 
scattering mechanism of metallic nanostructures similar to 
x-ray scattering by oscillating electrons. That is, free elec- 
trons on a non-planar conducting surface cannot move ho- 
mogeneously in response to the incident wave, thus forming 
inhomogeneous charge distributions. The oscillating charges 
then emit wavelets, and the interference between the wavelets 
may give rise to anomalous light scattering. The basic re- 
quirement here is free electrons, so this mechanism can ex- 
plain scattering from various conducting nanostructures, in- 
cluding perfect conductors and conductors with Re(e c ) > 
(see the TM' curve in Fig. 1 and the experiments and simula- 
tions of tungsten hole arrays with positive Re(e c ) in |U,[2(J). 
It is also applicable to non-periodic structures. For example, 
Figure 3 indicates that an isolated slit (or hole) can still be a 
single light source emitting divergent light. If the slit is sur- 
rounded by grooves (in which the cavity resonance may be 
different though), the grooves provide additional light sources 
that suppress the oblique wavelets (forming SSPs). Mean- 
while, wavelets emitted from the grooves along the back- 
ward direction offset specular reflection. Thus, transmission 
is enhanced. If grooves also exist on the exit surface, they 



again suppress oblique wavelets so that a collimated zero- 
order transmitted beam can be achieved [2]. As another exam- 
ple, it is known that when a nanowire is illuminated by a TM 
wave at one end, light can be "transferred" to the other end 
lf2ll [22tl . The common explanation is that light is transferred 
by CSPs on the wire, but the dispersion trend in l2lll is quite 
different from that of CSPs. Based on our mechanism, the in- 
cident wave causes inhomogeneous charges at the illuminated 
end that tend to propagate away due to the charge-netural ten- 
dency. This is similar to normal electricitiy transmission over 
metallic wires. The charge movement is then discontinued at 
the other end, leading to charge accumulation and oscillation 
that emit new light. The electrons can be bounced back, re- 
sulting in a Fabry-Perot-like charge pattern (that obviously do 
not need to obey the CSP principle). Here high conductivity 
can enhance the transfer efficiency, which is opposite to the 
CSP prediction that (nearly) perfect conductors do not sup- 
port CSPs. Overall, our mechanism reveals the fundamental 
of metallic nanostructure optics. 
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